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Abstract. In this study, a model of a Schottky-barrier carbon nanotube field- 
effect transistor (CNT-FET), with ferromagnetic contacts, has been developed. The 
emphasis is put on analysis of current-voltage characteristics as well as shot (and 
thermal) noise. The method is based on the tight-binding model and the non- 
equilibrium Green's function technique. The calculations show that, at room 
temperature, the shot noise of the CNT FET is Poissonian in the sub-threshold region, 
whereas in elevated gate and drain/source voltage regions the Fano factor gets strongly 
reduced. Moreover, transport properties strongly depend on relative magnetization 
orientations in the source and drain contacts. In particular, one observes quite a large 
tunnel magnetoresistance, whose absolute value may exceed 50%. 



PACS numbers: 85.75.Hh, 75.47.De, 42.50.Le, 73.63.Fg 
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1. Introduction 

In view of well-known size reduction problems of the conventional Si-based electronics, 
there have recently been intensive studies on new technologies based on nanostructured 
materials which are formed by organized growth and self-assembly methods pQ. A 
remarkable example of such self-organized structures are carbon nanotubes, which 
due to their fascinating physical properties have been studied for more than one 
and a half decade now and are believed to have commercial applications in the near 
future. In particular a nanotube-based transistor was first demonstrated in [2], and 
has since been studied extensively Carbon nanotubes are also very promising 

for spintronic applications, they may act as magneto-resistive switches [HI [7], and can 
maintain spin coherence over long distances. The latter is crucial for the tunnel (giant) 
magnetoresistance TMR (GMR) effect to be observable. In fact the TMR effect in CNTs 
was first reported in [8] , where the spin diffusion length of electrons flowing through a 
carbon tube was estimated to be on the order of routinely used contact separations 
in electric transport measurements. Subsequent studies confirm that the TMR effect 
in CNTs can be quite large, ranging typically from a few up to several tens percent 
([9]- [IS]). Here it is shown that also for the Schottky-barrier nanotube transistor the 
situation is likewise. 

Shot noise is another important phenomenon of present interest. It only appears 
in a non-equilibrium situation (finite source/drain voltage) and originates from the 
discreteness of the electron charge. In spite of there being a lot of papers on shot noise in 
nanostructures ([IB]- [IS]) (see [20] for a comprehensive review of recent advances), there 
have hitherto been no attempts to investigate shot noise in ferromagnetically contacted 
Schottky-barrier nanotube transistors (to the author's knowledge). 

The paper is organized in the following way: In Sec. 2 a short discussion of the 
adopted approach is presented, including the way the electrostatics has been dealt with, 
and basic equations concerning the non-equilibrium Grean's function method. Section 
3 is devoted to main results of the paper, i.e. current and shot noise dependences on 
the gate and drain/source voltages, as well as the impact of ferromagnetic electrodes on 
the above mentioned transport characteristics. Finally, Sec. 4 concludes the papers. 

2. Methodology 

A nanotransistor under consideration consists of a single wall carbon nanotube (CNT), 
described in terms of vr-orbital electrons, end-contacted to s-type itinerant-electron slabs 
of fee (111) crystallographical structure. The latter may be either paramagnetic or 
ferromagnetic. The system considered here is end-contacted, in contrast to other wide- 
spread geometries, like those of side-contacted or embedded ones. The device has been 
relaxed in order to find energetically favorable positions of interface atoms represented 
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by big and small spheres with diameters of = 2.5lA and dc = 1.42lA, for metal and 
carbon, respectively [TUl IH]. The studies are carried out within the framework of the 
tight-binding model and the non-equilibrium Green's function technique. The basic idea 
of this study is to adopt the methodology known for one- dimensional transistor models 
[211 [221 [23] . so as to make it useful for description of a Schottky-barrier nanotube 
transistor. This method bears much similarity to earlier methods [211 [251 126], except 
that here the calculations are performed in real space for CNT of finite length, and 
there is no need of using any ideal nanotube energy spectra (nor the Wentzel-Kramers- 
Brilloin, WBK, approximation). It should be however stated in this context that there 
is a great deal of theoretical studies on nanotube transistors (see e.g [27]- [31]), which 
successfully avoid reducing the Poisson equation to the 1-D problem, and handle thereby 
the electrostatics on a higher level than it is done here. Nevertheless, the present 
approach, while benefiting enormously from the mathematical simplicity, still leads to 
qualitatively correct results. The one- dimensional Poisson equation and its solutions 
are regarded here as spin-dependent (however spin indexes are skipped for brevity), and 
read [21] 

;Vix) + ^[Vg- Vix)] + ^—rpix) = 0, (1) 



dx^ A2 eoecNTA 



V{x) = / dx'v{x, x')p{x') + Ve,t{x), (2) 



^(x,x') = ^[e-l^-^-'l/^-e-(^+^')/^ + 

e"-^/^ (cosh ^ ~^ - cosh ) / sinh (3) 

V A A / A 



Vext{x) = [Vs sinh — h Vd sinh -]/ sinh - 

+ dx'v{x,x')VG, (4) 

where A = {R/2)^2{ecNTl^ox) ln(l + dox/R) + 1 is the effective screening length 
[23], ^CNT and R stand for the dielectric constant and the radius of the CNT, whereas 
the subscript ox refers to the coaxial gate oxide layer 5^02 {dox in thickness). The other 
symbols have the following meaning: Vs, Vd, Vg - source, drain and gate voltages; L - 
the length of the CNT; fs^o - Fermi function at Vs and Vd; A - the CNT cross-sectional 
area. 

As regards the non-equilibrium Green's function (GF) technique, a recursive 
method, similar to that of [32] was used, with obvious modifications consisting in 
including the potential profile V(x) (Eq. [2]) to the CNT Hamiltonian, and the Vs and 
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Vd voltages in the contacts. This procedure corresponds to the Hartree approximation, 
which is sufficient for the present room-temperature model, far beyond the Coulomb 
blockade and Kondo regimes. With this proviso one can write down equations for spin- 
projected current (I), zero- frequency noise power (S) and transmission matrix (T) as 
follows 



where trace (Tr) is taken over the orbital indexes, and = iijy^ — TJ^) (with 
a = S, D for source and drain, respectively). The first term of the noise power S is the 
thermal noise, it disappears at zero temperature, but dominates at finite temperatures 
and vanishing source/drain potentials. The Fano factor is defined here - assuming that 
both spin channels are independent (no spin relaxation processes) - in a usual way as 
F = {S^ + S^)/{2e{P + /^)), so that it is one for the Schottky shot noise (Poissonian 
limit). 

Instead of a rather common wide-band type approximation, the present method 
uses energy- dependent self-energies (S„(£')) which have been expressed in terms of: 
(i) recursively computed surface Green functions {ga{E)) of infinite fcc-(lll) contacts, 
and (a) the CNT/contact coupling matrices (vc), i.e T,a{E) = VcQaiE) vl- Here the 
retarded, G^, (advanced, G"^) GF is a matrix of rank equal to the number of atoms in 
the unit cell it refers to. The density matrix is defined as p = —i/{2ii) J dEG^{E) 
in terms of the lesser GF, which can be brought into the following form: G'^ = 
— /d(G"' — G"-) + i{fs — fD)G^TsG°'. While integrating G<, its first term contribution 
is found by integral over the contour in the complex energy plane, and the second one 
- by integrating along the real energy axis (using the Gaussian quadrature). The so- 
determined integral makes it possible to get an extra (excess or deficit) electric charge 
per CNT atom, with respect to the charge neutrality situation (no external potentials). 
The procedure is self-consistent because the potential V enters the diagonal of the 
Hamiltonian, which in turn, is present in the GF denominator. In order to make 
the recursive algorithm effectively work, and speed up the convergence, an additional 
assumption has been made, viz. the potential profile is only dependent on the unit cell 
number and remains constant within the given unit cell. 




dEUs-fu) Tr [T{E)l 



S = ^jdE {[fs{l - fs) + /d(1 - fo)] Tr[TiE)] 
+ {fs-fD)'Tr[T{E){l-T{E))]}, 



T = TsG'TdG", 



(5) 
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Figure 1. Spin-split surface densities of states (DOS) of a contact (thin line). In the 
antiparallel configuration the second contact has got interchanged \ sub-bands (no 
interchange in the parallel case). The thick curve is for the CNT DOS. Inset: Band 
bending at the interface for Vg > (short-dash line) and Vg < (solid line). 

3. Results and discussion 

Detailed computations have been carried out for a zigzag semiconducting CNT with 
a chiral vector (n, 0) and n=14. The radius R = an/(2-7r), and the energy gap 
A = 2^lt/{^/3n), where a=0.249 nm is the graphene lattice constant and t = 2.7 eV 
is the vr-orbital hopping integral. Additionally, it has been set: (^cnt — 1; ^ox — 4, 
dox = 2.5 nm. The CNT is I unit cells long (L = la\/3), and the calculations are 
performed for / = 100 (as well as for the ultra-short length / = 15). It is well known 
that the Schottky-barrier at a semiconductor/metal interface depends mainly on relative 
work functions (WF) of materials which form the junction [301 EH [33]. Here it is 
assumed that the metal WF is greater than that of the CNT, as e.g. for Pt, Au, 
Ni, Fe (as opposite to Ti), and locate the Fermi energies below the mid-gap. To be 
specific, the metal Fermi energy has been fixed at one-quarter of the band gap. On- 
and off-states of the FET are illustrated in Inset to Fig. [1], where the lower curves 
correspond to the maximum of the valence band and the upper curves correspond to 
the minimum of the conduction band. Accordingly, the present model-device represents 
an "Ep close to valence band maximum" case (cf. |1]). It is seen that, depending on 
the gate-voltage, the CNT-FET is either in the on-state (solid lines) or the off-state 
(short-dash lines). The former is due to the hole tunnelling through the Schottky- 
barrier, whereas the latter results from the fact that the transport energy-window falls 
within the energy gap. The location of Ep beyond the mid-gap of the isolated CNT 
introduces the electron-hole asymmetry, as shown recently [31] a similar effect may be 
caused by the presence of the inner uncontacted shell in a double-walled CNT, leading 
essentially to similar results for current-voltage characteristics (in the non-magnetic 
case). Figure [2] clearly shows that depending on boundary conditions imposed by the 
magnetic contacts [parallel alignment (P) vs. antiparallel alignment (^-P)] the current- 
gate voltage characteristics differ considerably from each other. Notably, the ON/OFF 
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Figure 2. The current vs. gate- voltage characteristics for the parallel alignment of 
the contact magnetizations (solid line) and the antiparallcl alignment (dashed line). 
The drain/source voltage is equal to 0.2V, length L=100 unit cells, and temperature 
T= 300°K. 



current ratio is also different for both the ahgnments. This is because the sub-threshold 
region reveals a conventional (JuUiere's type [35]) behavior with Ip clearly greater than 
Iapi whereas beyond this region both the currents differ far less from each other. Within 
the "non- JuUiere's" region the negative TMR apparently appears as a result of the shift 
of threshold Vq for electron conduction between P and AP configurations. Figure [3] 
presents the corresponding Fano factors. The shot noise of the CNT FET is Poissonian 
(F=l) in the sub-threshold region, and gets substantially reduced for elevated Vq and 
Vds- a still stronger reduction would be possible if the transport regime e.g. was either 
balhstic or of electron billiard type (with a short dwell time) [Ml [36]. As expected, 
magnetic conditions of the contacts, which strongly influence the current, have also a 
considerable effect on the noise. It is readily seen from Fig. [3] that the Fano factor 
corresponding to the P alignment, when compared to that of the AP alignment, is 
predominantly suppressed. The suppression takes place in a large region of the gate 
voltage, where current in the P configuration is greater than in the AP configuration (cf. 
Fig. [2]), in concord with the fact that F ~ 1//. It should be stressed that the transport 
regime considered here is phase coherent, although with noticeably reduced transmission 
due to the presence of Schottky-barriers. That is why the length effects, as shown in 
the Inset to Fig.3, are quite moderate, in contrast to what could be expected in the 
diffusive transport regime. Spin-dependent electron wave interferences are responsible 
for non-monotonic behavior of presented plots and accompanying sharp features. Similar 
features have been reported both in conductance and TMR [HI [131 EZ] as well as in 
the differential Fano factor [12]. The TMR is defined here in a "pessimistic" way as 
TMR = (Jp — Jap)//p, where Ip (Zap) is the current flowing through the system in the 
P (AP) magnetic configuration. The calculations are performed for spin polarizations 
in the ferromagnetic contacts equal to (n| — 'ni)/{ni + = 0.5 (ng. is the number 
of cr-spin electrons per atom). The results concerning TMR against Vg and Vds are 
presented in Figs. 3 and 4, respectively. The corresponding Insets clearly show that 
the TMR absolute value may be quite substantial and exceed several tens percent both 
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Figure 3. As in Fig. O but for the Fano factor. Note that reversal of the magnetic 
configuration from AP to P, for Vg < 0, may result in a considerable suppression of 
F. Inset: TMR for 2 different CNT lengths and Vds = 0.21/. 




Figure 4. Drain/source voltage dependence of I, TMR and F (main, left and right 
panels) for Vg = 0.3V, T=300K and L=15. In the interval ±2kBT around Vds = 
thermal noise dominates and / ^ 0, so F is ill-defined. 

in the sub-threshold region and at elevated voltages. Interestingly enough, TMR can 
assume large negative values (inverse TMR) and is strongly sensitive to both the gate- 
and Vd5- voltages. Another important observation is a striking asymmetry of the Fano 
factor and the GMR effect with respect to the sign of the voltage (p-channel vs. n- 
channel). These features are due to the aforementioned location of Ep with respect to 
the mid-gap, which makes the Schottky barrier for electrons higher than that for holes, 
and consequently the hole current (for negative voltages) is higher than the electron 
current (for positive voltages). It is noteworthy, that the oscillations present in Figs. [2]- 
S] are due to the size-dependent quantization of the CNT energy-levels, so the local 
maxima (minima) therein depend on to which extent the voltage-driven active channels 
and the transport energy-window fit together. This mechanism is also responsible for 
the occurrence of the inverse GMR, in fact it is similar to the mechanism based on 
resonant tunneling suggested in [12]. 
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4. Conclusions 



An approach has been developed to describe nanotube transistors, which combines 
a simplified (analytical) treatment of the electrostatics, and the state-of-the-art non- 
equilibrium GF technique. Unprecedented studies of electric transport through 
Schottky-barrier CNT transistor with ferromagnetic electrodes have been carried out. 
It has been shown that the tunnel ability of particles depends strongly on the alignment 
of contact magnetizations. Remarkably, the tunnel magnetoresistance can exceed 
50%, whereas the Fano factor is usually subjected to some extra suppression when 
the magnetization alignment gets reoriented from the antiparallel to the parallel 
one, improving thereby the noise-to-signal ratio. These two findings are relevant for 
prospective applications of nanotubes in spintronics. 
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